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ABSTRACT 

NGC 253 is a nearby spiral galaxy that is currently undergoing a nuclear burst of 
star formation in a 100 pc diameter region. We present spatially resolved 8-13 fj,m low- 
resolution spectra at four positions along the ridge of 8-13 fim emission. We find that 
the relative strengths of the ionic, dust emission features, and dust continuum emission 
vary with position in the galaxy but can be accounted for everywhere without recourse 
to extinction by silicates. The brightest mid-infrared peak (which is displaced from the 
nucleus) has elevated levels of both continuum and 11.1-12.9 /im "plateau" emission, 
indicative of dust heated within a photo-dissociation region. Spectra obtained over the 
course of 3 yr at the position of the brightest mid-infrared peak show no significant 
time variation. 
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1 INTRODUCTION 

NGC 253 is a nearby (D ~3 Mpc; Tully 1988) S(B)bc in- 
frared bright galaxy whose infrared luminosity (^2 xlO 10 
Lq; Telesco & Harper 1980) is powered by a ~100 pc di- 
ameter burst of star formation in the circumnuclear region. 
The infrared luminosity of NGC 253 exceeds that of the 
Milky Way by a factor of 2 with a larger relative fraction 
arising from the nuclear region. It thus provides an impor- 
tant proving ground to test our understanding of luminous 
infrared galaxies that are powered by compact starbursts. 
In these sources the higher luminosity per unit volume and 
the higher average density in the ISM may produce effects 
that are infrequently observed in star-forming regions in the 
Milky Way. 

The central regions of NGC 253 are strongly affected 
by dust extinction at visual wavelengths. It is now generally 
agreed that the nucleus of the galaxy lies near to the point 
marked TH2 (Turner & Ho 1985) in Fig. 1, some 3" (45 pc) 
away from the brightest infrared peak PI (Sams et al. 1994; 
Kalas & Wynn- Williams 1994; Pina et al. 1992; Keto et al. 
1993). The nature of PI is unclear. Based on its mid-infrared 
properties, Kalas & Wynn- Williams (1994) suggested that it 
is a very young supernova remnant embedded in dust, while 
Watson et al. (1996), who observed it with the Hubble Space 
Telescope, suggested that it contains a young star cluster. 

The 8-13 fim spectral region provides useful informa- 
tion on the conditions in starbursts through thermal dust 
continuum emission, dust emission features, and ionic emis- 
sion lines. The five strongest of the dust emission features 
occur at 3.3, 6.2, 7.7, 8.6, and 11.3 /im. The features are 



frequently attributed to polycyclic aromatic hydrocarbons 
(PAHs) (Puget & Leger 1989), but other laboratory ana- 
logues have also been proposed (e.g., Duley 1989; Sakata & 
Wada 1989; Ellis et al. 1994). PAH emission is widespread in 
the Galaxy. The features are observed in photo-dissociation 
regions (PDRs), reflection nebulae, and infrared cirrus, as 
well as in planetary nebulae. The ionic emission lines [Ar III] , 
[Siv], and [Nell], on the other hand, all arise in Hll re- 
gions and can be used as indicators of both the hardness 
and strength of the ionizing continuum radiation in obscured 
regions. 

Roche & Aitken (1985) used a Y.'6 aperture to obtain 
8-13 ixm spectra at two positions in the starburst in the nu- 
cleus of NGC 253. They found spectra that are broadly sim- 
ilar to those of many other star-forming galaxies, with con- 
tributions from dust continuum, PAH emission, and [Nell]. 
They discuss evidence for silicate absorption in their spec- 
tra, a topic that we re-examine in this paper. 

In this paper we present new 8-13 fim spectra at four 
positions in the starburst nucleus along the mid-IR ridge of 
emission with particular emphasis on the peak PI, for which 
we have obtained observations over a 3-years baseline. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

The data were obtained using the cooled grating 10 and 
20 spectrometer CGS3 (Cohen & Davies 1995) at the 
United Kingdom Infrared Telescope (UKIRT) on Mauna 
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Kea. We used low-resolution mode (-^ ~60) so that the en- 
tire 8-13 fj,m atmospheric window was observed. Two sepa- 
rate interleaved grating positions are required to fully sam- 
ple a spectrum with the 32-element linear detector array. 
Observations were performed by chopping with a 15" throw 
in the northwest-southeast direction at a frequency of 5 Hz 
for 10 s, followed by a beam switch of equal duration. No 
more than 10 such cycles were performed before observing 
at the second grating position. 

2.2 Pointing and astrometry 

Fig. 1 shows the positions of our spectroscopic observations 
labeled as P0, PI, P2a, and P2b. These observations were 
obtained on three nights, one each in 1992, 1993, and 1994. 
Table 1 is a log of the observations giving the date, aper- 
ture size in arcseconds for full width at 10 per cent power, 
per-detector element on-source integration times in seconds, 
and for P0, P2a, and P2b the relative offset from the PI ob- 
servation in arcseconds. In Fig. 1, only the 2'.'4 aperture is 
shown for PI. The central position of the CGS3 aperture 
with respect to the optical guide camera was established by 
using an automatic peak up routine on bright stars. In 1992, 
1993, and 1994 the stars were HR 188, HR 8892, and HR 
268, respectively, which are all within 10° of NGC 253. The 
telescope coordinate system was established by observing 
2 to 3 Carlsburg Meridian Circle stars near the galaxy. In 
1992, and 1994, we then slewed to PI, while in 1993, we 
peaked up manually on the PI using the integrated flux be- 
tween 8 and 13 /im. We estimate that the 1992 and 1994 
slew positions are accurate to l'/5, while the offsets from the 
position of PI are accurate to 0'.'5. 

The publication of the Hipparcos catalog has allowed a 
more precise determination of the position of the star used as 
an astrometric reference by Kalas & Wynn- Williams (1994). 
In Table 2 we give the proper-motion-corrected position for 
SAO 166579 based on the Hipparcos catalog. The position 
has not been corrected for parallax, which is the main source 
of the quoted error. The conversion to the FK4 system was 
performed using the routine BPRECESS provided with IDL 
version 5 using the epoch keyword. The position is consistent 
with the proper motion corrected position derived from the 
FK4 format PPM catalog (available from the NASA Astro- 
physics Data System), which is 0:02 east and 1'.'34 south of 
the SAO position used by Kalas & Wynn- Williams (1994). 
Table 2 also gives the revised positions for Kalas & Wynn- 
Williams's Peaks 1-3, the radio positions of TH2 and TH7 
(Turner & Ho 1985), and the positions of two red objects 
observed by Watson et al. (1996). In Table 2 we see that the 
revised position for Kalas & Wynn- Williams's Peak 1 is now 
consistent with TH7 and that their Peak 2 is now about 1 a 
distant from TH2. The positions of Kalas & Wynn- Williams' 
Peaks 1 and 3 are consistent with Watson et al.'s Bright blob 
and Spot o, respectively; if Peak 1 and their Bright blob are 
brought into coincidence, then the positions of Peak 3 and 
Spot a agree within (f.01 in a and 0'.'02 in S, within the 
estimated internal accuracy of the two observations. 

2.3 Data Reduction 

Our spectra have been flux calibrated against bright stars. 
In 1992, 1993, and 1994, we used HR 1457 (N=-3.03; Toku- 



naga 1984), HR 7001 (N=0.00), and HR 7949 (N=0.0; as- 
sumed), respectively. Wavelength-dependent variations in 
instrumental response and atmospheric transmission were 
removed by dividing the object spectra by stellar spectra 
and multiplying by a model Planck function. In 1992 and 
1993, all spectra were first corrected to 1 airmass using at- 
mospheric transmission curves based on data collected dur- 
ing each observing run. In 1994, the same star was used as 
a spectral and flux calibrator and was observed within 0.1 
airmass of the source. In 1992, HR 0617 was used for spec- 
troscopic calibration with an assumed temperature of 4780 
K; in 1993, HR 0188 was used with an assumed tempera- 
ture of 4690 K; and for the 1994 observations, the assumed 
temperature for HR 7949 was 4690 K. All the spectral stan- 
dards are spectral type K0 or earlier so that they should 
not be strongly affected by SiO absorption (Cohen et al. 
1992a). Deviation from our model Planck function due to 
wavelength-dependent H~ opacity is also likely to be small 
(Cohen et al. 1992b). We estimate our absolute photometric 
accuracy to be better than 10 per cent. 

Wavelength calibration was based upon observations of 
a Kr lamp as described by Hanner, Brooke & Tokunaga 
(1995). 

Due to year-to-year variations in instrumental setup, 
which resulted in slightly different wavelength scales for our 
spectra, we have resampled all our spectra to a common 
wavelength scale using Gaussian weighting with a = 0.06 
/im. Thus, the individual data presented below are nearly 
but not entirely independent of their spectral neighbors. We 
found no significant difference between the spectra obtained 
at PI over the 3 yr, so we have combined them to enhance 
the signal-to-noise ratio. Given the compact size < l'/2 (Keto 
et al. 1993) of PI between 8 and 13 /im, it is not surprising 
that spectra obtained with our two aperture sizes should 
agree, and this agreement can be taken as evidence that PI 
was reasonably well centered all 3 yr. 



3 RESULTS 

In Fig. 2 we present our spectra of the starburst region of 
NGC 253. The panels are labeled in accordance with Fig. 1 
and Table 1. We note that certain features are apparent in 
all the spectra. The brightest emission features are the [Ne II] 
forbidden line at 12.8 /im and the PAH emission feature at 
11.3 /im. The 8.6-/im PAH feature is also present in each 
spectrum. The continuum emission has a broad minimum 
near 10 /im; we shall argue in Section 4.1 that this 'contin- 
uum' is actually a combination of a true thermal continuum 
between 10 and 13 fim and the long wavelength shoulder of 
the 7.7-/im PAH emission feature. Each spectrum also shows 
signs of an elevated level of emission between 11.1 and 12.9 
/im. This so-called "plateau" emission has been previously 
observed in the Orion Bar photo-dissociation region (PDR) 
by Roche, Aitken & Smith (1989). There are no signs of the 
10.5-/im [Siv] or 9.0-/im [Arm] forbidden emission lines in 
any spectrum. 

As a check on our calibration, we have computed the av- 
erage flux density between 11.15 and 12.25 /im of the bright 
source PI. Our value of 4.7±0.1 Jy is in good agreement 
with the estimate of Pifia et al. (1992) of 4.8 Jy. 
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4 DISCUSSION 

Our discussion will first concentrate on two progressively 
more detailed fits to our spectra culminating with a table of 
measurements of spectral feature properties. We will then 
discuss the question of silicate absorption in NGC 253, fol- 
lowed by an analysis of the strengths of the ionic emission 
lines. Finally we discuss the possible origins of the contin- 
uum emission and the spatial variations seen in our spectra. 

We fit the features in our spectra in two stages. First 
we apply a two-component model to the whole of each spec- 
trum consisting of a theremal continuum underlying PAH 
emission. This procedure allows us to test whether or not 
silicate absorption is present in these spectra by examining 
the agreement of the observations with what the fits predict 
near 8 fim. Second, we decompose the emission above the 
thermal continuum longward of 11 fim into three compo- 
nents: PAH feature emission, [Nell] emission, and plateau 
emission. 



4.1 Initial fitting: Continuum and 11.3 (im feature 

Our initial procedure is to fit the 10-13 fim spectrum of 
each region to a combination of thermal blackbody and 11.3- 
fim PAH emission feature, and then compare this model's 
predicted 8-10 fim emission with the observations. The pro- 
cedure is similar to that employed for Arp 299C in our ear- 
lier paper (Dudley & Wynn- Williams 1993) except that the 
continuum is now estimated solely from the spectra them- 
selves rather than in conjunction with longer wavelength 
photometric data. For a PAH template we use the Orion 
Bar data of Roche et al. (1989) degraded in spectral resolu- 
tion to match the present data. The particular spectrum we 
use is that taken at Beckin et al.'s (1976) position 4, follow- 
ing Roche (1989) so that a prediction can be made of the 
strength of 8-10 fim emission. 

In Fig. 2, the solid lines are Planck function fits to the 
continuum between 10 and 11 /im and at 13.2 fim and the 
dotted line shows the Orion Bar spectrum scaled so that the 
intensity of Orion Bar emission around the 11.3-fim PAH 
feature between 10.9 and 11.7 fim matches that found above 
the continuum fits in our spectra. The shaded regions in 
Fig. 2 represent the reasonable range in the predicted 8-10 
fim emission based on the 10-13 fim fits. They are bounded 
on the bottom by the Orion Bar spectrum of Bregman et 
al. (1989) with the same scaling relative to the Roche et al. 
(1989) spectrum as that adopted by Roche (1989), and on 
the top by the same multiplied by a factor of 2. The shaded 
region reflects the observation that there is a factor of 2 
scatter in the ratio of the intensities of the 11.3 to 7.7+8.6 
fim features reported by Zavagno, Cox, & Baluteau (1992) in 
observations of galactic Hll regions. Of the sources they an- 
alyze, the Orion Bar is among those with the strongest 11.3- 
fim features relative to the 7.7+8. 6-/u,m features. As seen in 
Fig. 2, the simple model described in the last paragraph pro- 
vides an acceptable description of the gross features of all 
the spectra of NGC 253 we obtained. We therefore conclude, 
leaving aside the ionic [Nell] emission line, that the 8-13 
fim spectrum NGC 253 can be explained by a thermal con- 
tinuum rising towards longer wavelengths, plus a PAH spec- 
trum scaled to the intensity of the 11.3-/im feature. Thus 
the minimum near 9.7 /im need not be produced by silicate 



absorption; we discuss the limits on silicate absorption in 
Section 4.3. 

The color temperatures and 13.2-fim continuum flux 
densities of the four measured regions are given in rows 1-2 
of Table 3. As a check on the appropriateness of modeling 
the continuum in this manner, we extrapolate the results of 
the fits to 19.5 /im and compare the extrapolations with the 
data given by Pina et al. (1992). For IRS 1, they report a 

19.5 fim flux density of 24.8 Jy, which is within 6% of our ex- 
trapolation at PI of Fig. 1. For all the other positions in Fig. 

I, our extrapolations over predict the contours of their 19.5 
/jm deconvolved map. The largest over predictions (a fac- 
tor of 2) occur at P0 and P2b where we measure the lowest 
values of the color temperature. There are several possible 
explanations for this. First, notwithstanding Pina et al.'s 
efforts to deconvolve their map, there may remain impor- 
tant differences between their beam size and ours. Second, 
a blackbody may not be the best representation for the con- 
tinuum; a A -1 or A -1 ' 5 emissivity law gives a much better fit 
to their map. Third, some silicate absorption may have led 
to an underestimate of the continuum temperature at these 
positions. The silicate optical depth required (r a ii(9.7) ~ 0.7, 
with Tcont = 150 K) does not exceed the limits we set on 
the silicate optical depth in section 4.2. 

4.2 Detailed fitting: The plateau emission 

In the Orion Bar, the [Nell] emission, 11.3-pm PAH feature, 
and plateau all have different spatial distributions (Roche 
et al. 1989); we have chosen to fit these three components 
separately after first subtracting the thermal continuum. 

Our fitting method is quite similar to that employed 
by Roche et al. (1989) to analyse the Orion Bar. There are 
two main differences: (1) In the Orion Bar, a feature at 12.7 
/mi shows the same spatial variation as the 11. 3- /urn feature 
(Roche et al. 1989), but due to our lower spectral resolu- 
tion, the 12.7-/im feature cannot be independently measured 
here. We have therefore assumed a fixed ratio between the 
intensities of the 11.3- and 12.7-fim features of 2.4 based on 
the Orion Bar results. (2) We have chosen a plateau spectral 
profile that is stronger at longer wavelengths in Fa rather 
than the flatter profile employed by Roche et al. (1989). 
While our profile has the same spectral width as theirs, it 
rises linearly in wavelength and is twice as strong at 12.5 
fim than at 11.5 /im. This shape was chosen to match the 
shape of the plateau emission in the PI spectrum, which has 
the highest signal-to-noise ratio. 

We first fitted and subtracted the plateau by match- 
ing intensity between 11.6 and 12.2 /jm for the data and 
the profile. We then matched intensities and subtracted an 
assumed Gaussian profile centered at 11.3 fim and with a 
= 0.13 fim, where the matching occurred between 11.0 and 

II. 4 fim; a Gaussian centered on 12.7 fim with width a — 
0.26 fim was also subtracted as described previously. Finally 
the intensity of the [Nell] feature was measured between 

12.6 and 13.0 fim. These procedures left residuals that were 
not significantly different from zero except in the PI and 
the P2a spectra, where an apparent feature at 11.5 fim of 
intensities 1.3 ± 0.4 and 1.4 ± 0.4 xl0~ 15 W m~ 2 , respec- 
tively, measured between 11.3 and 11.7 fim remained. Since 
such a feature has been observed in some Galactic sources 
exhibiting PAH emission (Roche, Aitken & Smith 1991a), 
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these features could be real, but since we have modified the 
shape of our plateau model from that employed by Roche 
et al. (1989), we do not make any strong claim for this. 
This fitting procedure also produced a lower estimate of the 
11.3-/im intensity than the model used to produce Fig. 2. 
The largest reduction occurred for PI, where the new 11.3- 

intensity is about 80 per cent of that used in Fig. 2. 
About half of this reduction is due to plateau subtraction, 
and half is due to our effort to exclude flux from the appar- 
ent feature at 11.5 /mi. 

We summarize the results of the fits shown in Fig. 3 in 
rows 3-5 of Table 3. Rows 6 and 7 are limits on the intensities 
of two additional ionic emission lines: 8.99-/im [Arm], and 
10.5-/xm [Siv]. The equivalent widths of the 11.3-/tm PAH 
feature and the 12.8-/im [Ne II] line in rows 8 and 9 are based 
on the intensities in rows 3 and 5 and a model continuum 
that includes contributions from both the thermal contin- 
uum given in rows 1 and 2 and a plateau contribution, and 
in the case of [Nell], an assumed 12.7-/im feature contribu- 
tion. Row 10 gives the intensity ratios between the plateau 
and 11.3-/Ltm PAH features. The errors given in Table 3 are 
1 a while upper limits are 3 a. 

There are few spectra of galaxies available with 
sufficient signal-to-noise ratios to detect the 11.1-12.9 
/im plateau feature, so it is difficult to make comparisons 
with other galaxies, although plateau emission is clearly 
present in the two Infrared Space Observatory spectra of the 
antennae galaxy system presented by Vigroux et al. (1996). 

Most of what is known about the plateau emission is 
based on observations of the Orion Bar region. Roche et 
al. (1989) have made spatially resolved observations of the 
Orion Bar. Their observations of the ratio of the intensities 
of the plateau emission to 11.3-/itm feature varies between 
1.5 and 6.7. The lowest value is found on the visible Bar 
itself (Becklin et al.'s (1976) position 4), which also behaves 
to some degree as a reflection nebula, and the highest value 
is found at a position 20" south, well within the dense part 
of the molecular cloud. At the position where the plateau 
emission is the strongest (5" south of position 4), the ratio 
has a value of 2.7. The range of values for this ratio listed in 
Table 3 for NGC 253 falls within the range of values we have 
just calculated. Of particular interest are the high values of 
the ratio seen at PI and P2a that can only be approximated 
by the southern most values in the Orion Bar. These values 
may show that conditions similar to those occurring in the 
deep layer of the Orion Bar PDR are more prevalent than 
those occurring adjacent to the ionization front. 



4.3 Limits on Silicate Absorption 

Aitken & Roche (1984) have proposed a method of deter- 
mining the amount of silicate absorption required for spectra 
that include PAH emission based on observations of galac- 
tic planetary nebulae and employing the equivalent width 
of either the 11.3- or 8.6-/im features. Given the factor of 
2 dispersion in the observed 11.3- to 7.7+8.6-/im features 
intensity ratio in galactic Hll regions, we prefer to set an 
upper limit on the amount of true thermal continuum that 
might be present at 8 /im and thus an upper limit on rsn(9.7) 



where 

(Aitken & Jones 1973). We therefore take the difference be- 
tween the lower boundary of the shaded region at 8 /im (less 
our model continuum at 8 /im) and the spectra as a rough 
limit on the amount of continuum that can contribute to the 
flux at 8 /jm. This value, along with the model continuum 
flux at 9.7 and 13 pm, combine to give an upper limit on 
7"Sii(9.7) in the equation of Aitken & Jones (1973). These 
limits are given in row 11 of Table 3. 

We find that our upper limits on rsn(9.7) are in agree- 
ment with the tsu(9.7) reported by Roche & Aitken (1985) 
of 1.4 ± 0.3 based on the method of Aitken & Roche (1984) 
for NGC 253, given the very different aperture sizes we have 
employed. We note that these limits on rsn(9.7) apply to sil- 
icate absorption by cold dust, but not to the optical depth in 
dust that is emitting between 8-13 /im. The limit for PI is 
consistent with the amount of extinction estimated by Wat- 
son et al. (1996) based on V — I colors under the assumption 
that A v = 15 x r Si i(9.7) (Aitken 1981). 

4.4 Ionic Lines 

The presence and absence of ionic lines in our spectra allow 
us to derive properties of the ionized gas within NGC 253. 

In all four of the spectra shown in Fig. 2, the [Ne II] emis- 
sion line is clearly detected, but the [Ar in] and [S iv] lines 
are not. Of these two lines the limits on I([Siv]) in Table 3 
are probably more reliable given the systematic errors intro- 
duced by both the presence of the shoulder of the 7.7+8.6- 
/xm features and the adjacent telluric O3 absorption that 
could affect the continuum estimate for the 8.99-/im [Arm]. 
We therefore concentrate on what can be learned from the 
limits on I([Siv]) emission. The ratio of I([Siv]) to I([NeIl]) 
is sensitive to the hardness of the Ly cont radiation field ion- 
izing the H II regions producing the emission and thus to the 
presence of very hot stars. Assuming a number abundance 
ratio of S to Ne of 0.27, n e ~ 10 3 -10 4 cm" 3 , and Ly cont 
fluxes following the models of Kurucz (1979), Rubin (1985) 
gives ratios of I([Slv])/I([Nell]) of 0.02 and 1 for exciting 
stars of T e g — 33,000 and 36,000 K, respectively. Our up- 
per limits on I([Siv]) imply that I([S Iv])/I([Ne II]) < 0.33, 
0.09, 0.14, and 0.13 for each of the spectra P0, PI, P2a, and 
P2b, respectively. Thus the nondetection of the [Siv] line 
suggests that the bulk of the stars producing ionizing radia- 
tion have Tog < 36,000 K. The strongest limit is seen at PI, 
where a finer grid of models than those examined here would 
be useful. The limit on T c ff holds even if rsn ~ 1.2, which 
would increase our upper limits on I([S Iv])/I([Ne II]) by a 
factor of 1.7 based on the Mathis (1990) extinction curve 
with R ~ 3 supplemented by a p Cep silicate profile match- 
ing our spectral resolution (see Dudley & Wynn- Williams 
1997 for details.) 

The [Ne II] emission may be compared with observations 
at other wavelengths. Puxley & Brand (1995) measure a Br7 
flux of 0.18 x 10~ 15 W m~ 2 at PI using a 3" slit oriented 
along the galaxy's major axis. This gives I(Br7)/I([NeIl]) 
= 0.015 which is within 25 per cent of the prediction of 
Dudley & Wynn- Williams (1993) of 0.011 but is quite low 
compared to the models of Rubin (1985) in the range of T e s 
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suggested by the upper limit on I([Siv]), which predict a 
ratio of 0.048. Different assumptions about Ne abundance, 
Brackett to Balmer series intensity ratios, and the [Ne II] col- 
lision strength can account for about half of the discrepancy 
between the two predictions. On the other hand, applying 
the reddening deduced by Watson et al. (1996) towards PI 
to the prediction of Rubin (1985) gives good agreement with 
the observed ratio. Applying the I([NeIl]) to I(Br7) conver- 
sion factor of Dudley & Wynn- Williams (1993) followed by 
a further conversion using the Br7 intensity to 6-cm free- 
free flux density ratio of 8 x 10 11 Hz from Wynn- Williams 
(1982) to the observed [Nell] intensity predicts a 6-cm flux 
density of 17 mjy, which is not too different from the ob- 
served value of about 10 mjy for TH 7 (Ulvestad & An- 
tonucci 1997). The 6 to 2 cm spectral index of this source is 
a — —0.24 (S v oc f a ), intermediate between optically thin 
free-free and synchrotron indices, suggesting that both pro- 
cesses contribute to the radio emission. 

For P2a, at the nucleus of the galaxy, the predicted 
Br7 flux based on the [Ne II] flux (Dudley & Wynn- Williams 
1993) is O.lOx 10~ 15 W m~ 2 , which is a factor of 2 lower than 
the 0.24x 10~ 15 W m~ 2 observed by Puxley & Brand (1995). 
This may indicate that there is an ionizing source associated 
with the nonthermal radio source (TH2) that produces Br7 
emission that is not mirrored in the [Ne II] emission, possibly 
because the ionized gas density is too high (n e > 10 5 cm~ 3 ) 
to produce as much [Nell] emission. This speculation is not 
supported, however, when using the ratio of Rubin (1985) 
which overpredicts the observed Br7 flux by a factor of 2. As 
in the case of PI, this effect might also be accounted for by 
reddening. Both comparisons are beset by possible aperture 
differences, an unknown Ne abundance, and uncertainties in 
determining the [Nell] intensity due to our subtraction of 
an assumed 12.7-/im PAH feature so that only firm conclu- 
sion that can be drawn is that the Br7 and [Nell] are in 
agreement within a factor of 4 under the assumption that 
they both arise from H II regions. 

The 6-cm radio continuum in P2a is about 40 mjy (Ul- 
vestad & Antonucci 1997), stronger than the free- free emis- 
sion predicted from the intensities of either [Nell] or Br7, 
which amount to about 13 mjy. This is consistent with the 
statement of Turner & Ho (1985) that the 2-cm radio con- 
tinuum emission from TH2 must be due to nonthermal emis- 
sion since its brightness temperature of 10 5 K is too high to 
be due to free-free emission from Hll regions. On the other 
hand, the amount of free-free emission that we would pre- 
dict by substituting the [Ne II] to Br7 conversion factor given 
by Rubin (1985) but retaining the Br7 to 6-cm flux den- 
sity factor adopted from Wynn- Williams (1982) gives good 
agreement with the observed 6-cm emission. 

4.5 Continuum emission 

In standard Case B recombination theory, radio free-free 
emission, Brackett series recombination line emission, and 
[Nell] emission all trace the number of ionizing photons 
with little sensitivity to electron temperature and density 
for the values typical of H II regions. Thus, [Ne II] emission 
may be used to predict the amount of grain emission due to 
trapped Lya heating in H II regions while avoiding the pos- 
sible problems of extinction at short wavelengths and non- 
thermal emission at radio wavelengths. Taking the Br7-to- 



[Ne II] intensity ratio predicted by Dudley & Wynn- Williams 
(1993) and a Br7 intensity to 6-cm free- free flux density ra- 
tio from Wynn- Williams (1982), we predict a free-free flux 
density at 6 cm of between 5 and 17 mjy for the four spec- 
tra. Using the continuum temperatures given in Table 3, wc 
may then use formula A4 of Genzel et al. (1982) to predict 
the expected flux density at 13.2 /jm due to trapped Lya 
heating of dust and compare these predictions with the ob- 
served values given in Table 3. We find that the observed 
continuum exceeds these predictions by large factors of 30- 
50. We therefore conclude that Lya heating fails to account 
for the 10-13 fim continuum in any of our spectra. 

4.6 Spatial Variations 

Before discussing the possible reasons for spatial variations 
among our spectra of NGC 253 we need to review the status 
of the compact source PI. 

In this paper we have already shown that the 8-13 
//m spectrum of PI is consistent with a combination of ther- 
mal dust emission and Orion Bar-like emission along with 
[Ne II] emission similar to that seen in starburst galaxies (see 
Fig. 2). Further, we observed no time variation in the 8-13 
/mi spectrum over the 3 yr that we observed it. In view of 
the absence of time variations, and of Watson et al.'s (1996) 
HST detection of a star cluster spatially coincident with PI, 
the suggestion of Kalas & Wynn- Williams (1994) that PI is 
a dust-embedded recent supernove is probably incorrect. 

Our spectrum of PI can provide a new estimate of the 
1-1000 /im luminosity of the source based on the color tem- 
perature and flux of the continuum. A 145-K blackbody with 
a 13.2-pm flux density of 6.7 Jy at a distance of 3 Mpc has a 
luminosity of 1.6x 10 9 L Q . We recall the agreement noted in 
section 4.1 between this extrapolation and the 19.5 /im data 
of Pina et al. (1992). This luminosity estimate for PI is in 
good agreement with the estimate of Watson et al. (1996) 
of 1.8xl0 9 Lq which is based on their estimate of the frac- 
tion of the flux due to PI in the 12.5-pm map of Keto et 
al. (1993) applied to the Lio-350//m for the inner 30" of 
NGC 253 given by Telesco & Harper (1980). However, our 
estimate of the flux due to PI in the Keto et al. (1993) map 
is a factor of 2.3 larger than that of Watson et al. (1996), 
so the agreement is largely fortuitous. By either estimate, 
PI provided ~8 per cent of the infrared luminosity of the 
galaxy. 

The 11.3-/iim feature is clearly present in all our spectra 
but shows a large variation in EW. This is consistent with 
the observations of the related 3.3-/mi feature reported by 
Kalas & Wynn- Williams (1994), where it was found to be 
relatively weak compared to the continuum at PI. The EW 
of the 11.3-/xm feature ranges over a factor of 4 and is anti- 
correlated with the continuum temperature, while the ratio 
of the intensities of [Nell] to 11.3-/im feature cover a factor 
of 2 and all fall within the range reported by Roche et al. 
(1991b) for 19 galaxies showing both features. 

The young star cluster model for PI leads, at least qual- 
itatively, to a plausible explanation for the differences seen 
between the spectum at PI (and to some extent P2a) and 
elsewhere in the galaxy. Because the interstellar radiation 
field is presumably stronger within the compact PI star 
cluster than the interstellar radiation field elsewhere, reg- 
ular dust grains will be heated to higher temperatures as is 
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observed (Table 3). Since we have already shown that Lya 
heating is insufficient to heat the dust, we need to invoke di- 
rect dust heating by starlight. The prominence of the plateau 
emission in the spectrum of PI suggests that emission that 
arises in PDRs is particularly important in PI, so we sus- 
pect that the dust that gives rise to the continuum emission 
is likely to come mainly from the outer layers of molecu- 
lar clouds exposed to the starlight. The relative importance 
of large-grain (equilibrium) heating and small-grain (tran- 
sient) heating is difficult to establish solely with the new 
data presented here. 
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5 CONCLUSIONS 

We have presented spatially resolved 8-13 /im spectroscopy 
of the nuclear region of the starburst galaxy NGC 253. We 
find that high-mass star formation explains the observed 
spectra through a combination of PAH emission, [Ne n] emis- 
sion, and thermal continuum with little need to invoke ab- 
sorption by cold dust greater than Ay ~15. 

Based on the revised astrometry of Kalas & Wynn- 
Williams (1994) we identify PI with an optical red object 
observed by Watson et al. (1996) and with the 2-cm radio 
source TH7 observed by Turner & Ho (1995). 

We report the measurement of a low-contrast but 
high- intensity plateau of PAH emission in the 11.1-12.9 
/im range, and identify it with emission observed in PDRs 
in the vicinity galactic Hll regions. 

We set a limit on the effective temperature of the in- 
tegrated emission from stars sampled by our spectra to be 
less than 36,000 K based on the comparison of I([Nell]) and 
upper limits on I([S IV]) under the assumption that the dust 
optical depth is not very large at 10.5 /jm. 

We compare I([Nell]) with I(Br7) and radio continuum 
emission from PI, and find reasonable agreement with stan- 
dard recombination theory, while for P2a the radio contin- 
uum emission is probably too strong to be due to Hll re- 
gions. 

We examine the relative intensities of 11.3-/im PAH 
emission feature, plateau emission, and the [Nell] emission 
line, and the continuum intensity and color temperature and 
find that while the [Nell] and 11.3 /im PAH emission are 
roughly proportional, as observed in other starburst galax- 
ies, the 11.3 /im PAH feature and plateau emission show 
a larger range of variation with respect to each other. We 
suggest that enhanced plateau emission is due to the pre- 
dominance of conditions more similar to the depths of PDRs 
than regions adjacent to ionization fronts. 
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Table 1. Observing Log for NGC 253 



Date Aperture 


Sources 


Time 


Offset 




(UT) (" ) 


Observed 


on Source (s) 


from PI 




1992 Oct 15 3.26 


PI 


zUU 








P2b 


500 


2'.'7E 27 ON 




1993 Aug 31 2.40 


P0 


500 


l'.'8W 2'.'0S 






PI 


350 








P2a 


500 


1"8E 2V ON 




1994 Aug 28 3.26 


PI 


100 






Table 2. Positions of Bright Peaks in 


NGC 253 






Publication 


A 


Source Name 


Position (1950) 








a = 00 h 45 m 


<5 = -25° 33' 


Kalas & Wynn- Williams 


4.8 ^m 


Peak 1 


05=632 ± 0.037 


41'.' 34 ± 0.5 


Turner & Ho 


2 cm 


7 


05=626 


41'.' 23 


Watson et al. 


0.6 fim 


Bright blob 


05=58 ± 0.02 


41" 08 ± 0.3 


Kalas & Wynn- Williams 


4.8 ^m 


Peak 9 


05=782 ± 0.037 


39'.' 50 ± 0.5 


Turner & Ho 


2 cm 


2 


05 s 795 


38'.' 95 


Kalas & Wynn- Williams 


4.8 lira 


Peak 3 


06=08 ± 0.037 


35'.' 89 ±0.5 


Watson et al. 


0.6 lira 


Spot a 


06=03 ± 0.02 


35'.'61 ± 0.3 


Hipparcos 




SAO 166579 


a = 00 h 45 m 


8 = -25° 40' 








7=4638 ± 0.0004 


3" 677 ±0.005 



Table 3. Spectral Measurements of NGC 253 

Position 



Property 


P0 


PI 


P2a 


P2b 


T cont [K] 


119l 9 8 


145t 2 


13611 


1191^ 


F„(13.2 /ira) [Jy] 
Iii.3 M m [lO" 15 W m- 2 ] 


1.7 ± 0.1 


6.7 ± 0.1 


4.3 ± 0.1 


3.3 ± 0.1 


3.3 ± 0.4 


6.0 ± 0.3 


6.7 ± 0.3 


5.9 ± 0.4 


I[Ncii] [HT 15 W m- 2 ] 


3.9 ± 0.8 


12.2 ± 0.8 


9.0 ± 0.8 


10.9 ± 0.9 


Ill.l-12.9 M m [10- 15 Wm- 2 ] 


6 ± 2 


30 ± 2 


25 ± 2 


17 ± 2 


I[Ariii] [10~ 15 W m- 2 ] 


<1.3 


<1.1 


<1.3 


<1.3 


I [S iv] [10- 15 W m- 2 ] 


<1.3 


<1.1 


<1.3 


<1.3 


EW 11.3 H 


0.22 ± 0.02 


0.075 ± 0.004 


0.14 ± 0.01 


0.19 ± 0.01 


EW [Neil] [urn] 


0.13 ± 0.03 


0.09 ± 0.01 


0.11 ± 0.01 


0.17 ± 0.01 


Ill.3Mm 


2.0 ± 0.6 


5.0 ± 0.5 


3.8 ± 0.3 


2.9 ± 0.3 


r Si i(9.7) 


<1.3 


<0.7 


<0.8 


<1.6 
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Figure 1. The positions and apertures for the spectra presented 
in this paper are shown as open circles and labeled P0, PI, P2a, 
and P2b. The filled circles are the positions of compact radio 
sources measured by Turner Ho (1985), with the larger two la- 
beled according to their table 2. Kalas & Wynn- Williams's Peak 
3 is also labeled. This map is adapted from the 4.8 (im map of 
Kalas & Wynn-Williams (1994). 



Figure 2. The 8—13 (im spectra at four positions in the starburst 
nucleus of NGC 253 are labeled as in Table 1. Large filled circles 
represent observed data points, while small filled circles, when 
present, indicate boxcar smoothed data, where the smoothing ra- 
dius is 0.1 fim. The solid lines are blackbody fits between 10 and 
11 (im and at 13 ^m. The temperature of the fitted blackbody is 
indicated in each panel. The dashed lines show a fit of the Orion 
Bar spectrum (Roche 1989) between 10.9 and 11.7 fim, while the 
shaded regions indicate the expected range in the strength of the 
7.7 fim feature based on the strength of the 11.3 fim feature. 



Figure 3. The continuum subtracted 10—13.3 /im spectra cor- 
responding to Fig. 2 are presented. The filled circles are the 
continuum-subtracted data displayed on linear scales in Fx and 
wavelength, with 1 a errors that do not include uncertainties in 
the continuum subtraction. The solid lines give the total model 
fits for each spectrum. The dot-dashed line in the P0 spectrum 
indicates the contribution of the plateau and 12.7 /im feature that 
are subtracted from the [Ne II] flux. 
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